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Abstract. A review of meson emission in heavy ion collisions at incident 
energies from SIS up to collider energies is presented. A statistical model 
assuming chemical equilibrium and local strangeness conservation (i.e. strangeness 
conservation per collision) explains most of the observed features. 

Emphasis is put onto the study of K + and K~ emission at low incident 
energies. In the framework of this statistical model it is shown that the 
experimentally observed equality of K + and K~ rates at "threshold-corrected" 
energies y/s — y/sth i s due to a crossing of two excitation functions. Furthermore, 
the independence of the K + /K~ ratio on the number of participating nucleons 
observed between SIS and RHIC is consistent with this model. 

It is demonstrated that the K~ production at SIS energies occurs 
predominantly via strangeness exchange and that this channel is approaching 
chemical equilibrium. The observed maximum in the ii"+/7r+ excitation function 
is also seen in the ratio of strange to non-strange particle production. The 
appearance of this maximum around 30 A-GeV is due to the energy dependence 
of the chemical freeze-out parameters T and /i g . 



1. Introduction 

Central heavy ion collisions at relativistic incident energies represent an ideal tool to 
study nuclear matter at high temperatures. Particle production is - at all incident 
energies - a key quantity to extract information on the properties of nuclear matter 
under these extreme conditions. Particles carrying strangeness have turned out to be 
very valuable messengers. 

A specific purpose of this paper is the presentation of the evolution of strange 
particle production over a large range of incident energies. Many results are 
shown together with a theoretical interpretation. The attempts to describe particle 
production yields with statistical models [1-8] have turned out to be very successful 
over this large domain of incident energies. 



2. General Trends 



2.1. Production of pions and kaons from SIS to RHIC 

At incident energies around 1 A-GeV pion and kaon production is very different: Pions 
can be produced by direct NN collisions in contrast to kaons. The threshold for K + 
production in AW collisions is 1.58 GeV and only collective effects can accumulate the 
energy needed to produce a K + together with a A due to strangeness conservation. 
The threshold for K~ production is even higher (2.5 GeV) as they are produced as 
K + K~ pairs. These conditions lead to very different yields for the various mesons 
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Figure 1. Spectra of positively and 
negatively charged pions and kaons 
measured in central collisions of Au+Au 
at 1.5 A-GeV. Preliminary results. 
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Figure 2. Spectra of positively charged 
pions in the center-of-mass frame close 
to midrapidity for Au+Au (upper part) 
and C+C (lower part). A deviation from 
a Boltzmann shape is seen at all incident 
energies and for both collision systems. 



as demonstrated in Fig. |l| showing some spectra from central Au+Au collisions at 1.5 
A-GeV. The yields of pions are much higher than for kaons; positively and negatively 
charged pions differ due to the N/Z ratio of Au. The yield of K + is by nearly two 
orders of magnitude higher than the one of K~ due to the different thresholds. At 
higher incident energies all these differences diminish as will be discussed. 

Pion spectra of excellent quality are now available or will appear soon [9- 
10]. Figure || shows as an example spectra of positively charged pions as d 3 a/dp 3 
(Boltzmann representation) for C+C and Au+Au collisions at different incident 
energies. The spectra are measured close to mid rapidity. 



Survey of experimental data 



3 




Figure 3. and K~ spectra measured at midrapidity for Au+Au (upper part) 
and C+C (lower part) at various incident energies |l -'[ |l3[ . The spectra of K~ 
from Au+Au collisions are preliminary. 



All spectra exhibitconcave shapes in this representation deviating from a 
Boltzmann distribution which would be a straight line. Yet, even in a thermal 
condition one is not expecting a pure Boltzmann distribution as pions originate both 
from "free" pions and from resonance decay after freeze out. This type of shape 
is observed up to the highest incident energies and the above interpretation holds 
qualitatively. For a quantitative description with these two components still some 
work has to be done. 

Spectra of K + and K~ from mass-symmetric systems C+C and Au+Au at 
incident energies from 0.6 to 2.0 A-GeV measured at midrapidity are shown in Fig. [| 
They exhibit Boltzmann shapes in all cases ]l2| , [l3| . Their inverse slope parameters 
increase monotoneously with incident energy and the heavier system exhibits harder 
spectra than the light system at the same incident energy. 

The interaction of K + and K~ with nuclear matter is very different: Due to their 
s content K + cannot be absorbed, while K~ can easily be absorbed on a nucleoli 
converting it into a A. This difference makes the K + to be messengers of the early 
stage of the collision. Therefore K + are ideal probes for this dense stage and allow to 
extract the stiffness of the nuclear equation of state [|l3| . The basis of these studies is 
the ratio of K + measured in Au+Au and in C+C collisions as shown in Fig. [l|. This 
subject is presented in detail in the talk by C. Sturm. 

Figure |3| demonstrates that the yield of K + is much higher than the one of K~ 
as mentioned already. This is caused by the different production thresholds for K + 
and K~ . At AGS energies the ratio of K + / K~ has decreased to about 5 JlJ] and is 



as low as 1.16 at RHIC [p_5|. This trend is summarized in Fig. 



5. The dashed line 



represents the results of calculations using a statistical model [|f| 5 
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Figure 4. Upper part: Multiplic- 
ities of K + and pions per A as a 
function of incident energy. The 
pion data include charged and 
neutral pions. The lines represent 
fits to the data (see Jl3||), Lower 
part: Ratio of the multiplicities 
per A (Au+Au over C+C) as a 
function of incident energy. 
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Figure 5. Ratio of positively and 
negatively charged kaons for various 
collisions systems from SIS up to RHIC. 
Preliminary results. The dashed line 
represents the results of calculations 
using the statistical model. 



2.2. Interpretation within a statistical model 

Pions and K + exhibit a further very pronounced contrast: While the pion multiplicity 
per number of participating nucleons A part remains constant with A part , the K + 
multiplicity per A part rises strongly (Fig. ||). The latter observation seems to be 
in conflict with a thermal interpretation, which - in a naive view - should give 
multiplicities per mass number A being constant. 

Usually, the particle number densities or the multiplicities per A part , here for 
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pions, are described in a simplified way by a Boltzmann factor 



A 



exp 



part 



T 



with the temperature T and the total energy < E n >. 

The production of strange particles has to fulfil strangeness conservation. The 
attempt to describe the measured particle ratios including strange hadrons at AGS 
and SPS using a strangeness chemical potential is quite successful jj], g, ff[ |j. 
However, this grand-canonical treatement is not sufficient if the number of produced 
strange particles is small. Then a statistical model has to take care of local strangeness 
conservation in each reaction as introduced in Jl6[ | . This is done by taking into account 
that e.g. together with each K + a A or another strange particle is produced: 
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where T is the temperature, fj,B the baryo-chemical potential, gi the degeneracy 
factors, V the production volume for making the associate pair (see j|, ||) and -E; 
the total energies. We note that this volume is not identical to the volume of the 
system at freeze out. The volume parameter V is taken as ryA part with a common 
ry = 1.07 fm for all systems and all incident energies. 

This formula, simplified for demonstration purposes, neglects other combinations 
leading to the production of K + as well as the use of Bose-Fermi distributions, which 
are all included in the computation. The corresponding formula for K~ production 
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is similar, but does not depend on ^tg. This point will become important later on. 

These formulae lead to a reduction of K + and if ~ yields as compared to the 
numbers calculated without exact strangeness conservation [Q, |5j. Two extreme 
conditions can be seen from these equations. In the limit of a small number of 
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strange particles the additional term (due to the parameter V) leads to a linear rise of 
M K + / A part , while M^/A part remains constant. This is in very good agreement with 
the experimental observations shown in Fig. ^|. For very high temperatures or very 
large volumina, the terms in brackets approach unity (see Ref. resulting in the 
grand-canonical formulation. This is much better seen in the exact formulae using 
modified Bessel functions |4|, |j| . 

At low incident energies, the particle ratios (except 77/710) are well described using 
this canonical approach Q. Surprisingly, even the measured K + /K~ ratio is described 
and this ratio does not depend on the choice of the volume term V. It should be noted 
that the statistical model uses nominal masses of the particles while some transport 
calculations have to reduce the K~ mass (as expected for kaon in the nuclear 
medium) in order to describe the measured yields. 

Before comparing the data with the calculations in detail, a summary of the 
measurements by the KaoS Collaboration is given. These results have attracted 
considerable interest as in heavy ion collisions the K~ yield compared to the K + 
cross section is much higher than expected from NN collisions [|l2| |20| . This is 
especially evident if the kaon multiplicities are plotted as a function of ^/s — ^Jsth 
where y/sth is the energy needed to produce the respective particle in NN collisions 
taking into account the mass of the associately produced partner. To produce a K + in 
NN collisions ^/Jth = 2.548 GeV and a K~ y/s^h = 2.87 GeV. The obvious contrast 
between NN and AA collisions, shown in Fig. ^, has lead to the interpretation of the 
results by in-medium properties which cause e.g. a lower threshold for K~ production 
when produced in dense matter ]l9| . The observed difference between NN and AA 
collisions alone is not sufficient to conclude on properties of kaons in matter. In heavy 
ion collisions, kaons can be produced by other channels, e.g. ttA — » K~ N which are 
not available in NN collisions. Only by using detailed transport-model calculations 
one might conclude on new properties of kaons in matter J19[ . 
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Figure 7. Measured K + and 
K~ yields in heavy ion and NN 
collisions as a function of y/s — 
^JE H @ < A par t > is 

A/2 forneavy ion data and 2 for 
NN collisions. 



It is therefore of interest to see how the results of the statistical model appear 
in a representation where the K + and K~ multiplicities are given as a function of 
\fs ~ -\/sth- Figure || demonstrates that at values of yfs — ^/Sth less than zero the 
excitation functions for K + and K~ cross leading to the observed equality of K + 
and K~ at SIS energies. The yields differ at AGS energies by a factor of five. The 
difference in the rise of the two excitation functions can be understood by the formulae 
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given above. The one for K + production contains — fig) while the other has E K + 
in the exponent of the second term. As these two values are different, the excitation 
functions, i.e. the variation with T, exhibit a different rise. 
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Figure 8. Calculated K + /A par t and 
K~ /Ap ar t ratios in the statistical model 
as a function of ^/s — ^fsth for Ni+Ni 
collisions. The points are resu 
Ni+Ni collisions at SIS energies 
and Au+Au at 10.2 A GeV (AGS} II 1 
At AGS energies the influence of th 
system mass is negligible. 




Furthermore, the two formulae predict that the K + /K ratio for a given collision 
should not vary with centrality as V cancels in the ratio. This has indeed been observed 
in Au+Au/Pb+Pb collisions between 1.5 A-GcV and RHIC energies |s| [H| ||, [lj, || 
as shown in Fig. This independence of centrality is most astonishing as one 
expects at low incident energies an influence of the different thresholds and the density 
variation with centrality. For instance at 1.93 AGeV the K + production is above and 
the K~ production below their respective NN thresholds. 
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Figure 9. The K+/K~ ratio 
appears to be constant of a function 
of centrality from SIS up to RHIC 
energies. The dotted lines represent 
the predictions of the statistical model. 
Data from |L4| 



Transport-model calculations show clearly that strangeness equilibration requires 
a time interval of 40 - 80 fm/c |23|, On the other hand statistical models assuming 
chemical equilibration are quite successful in describing the particle yields including 
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strange particles. 

In the case of K + production, no strong absorptive channel seems to be available 
which could lead to chemical equlibration. For K~ production the situation is quite 
different. At low incident energies strange quarks are found only in a few hadrons. The 
s quark is essentially only in K + , while the s quark will be shared between K~ and 
A (or other hyperons). This sharing of the s quark might be in chemical equlibrium 
as the reactions 

7T° + A ^ p + K~ or 7r~ + A ^ n + K~ 

are strong and have only slightly negative Q-values of -176 MeV. 

The idea that the K~ yield is dominated by strangeness exchange via the 7r~ + A 
channel has been suggested by p5| and has been demonstrated quantitavely in a recent 
theoretical study p6[| . The direct K + K~ pair production via baryon-baryon collisions 
has negligible influence as these K~ are absorbed entirely. In these transport-model 
calculations the strangeness exchange is approaching equilibrium but does not fully 
reach it p6fl . For details see the talk given by C. Hartnack. 

If these reactions are the dominating channels, they might reach chemical 
equilibration, e.g. the rates for producing K~ and for absorbing K~ arc equal. Then 
the law of mass action can be applied giving for the respective concentrations p7[ 

K-] • [A] _ „ 



[K~] ■ N 



As the number of K relative to A is small, [A] can be approximated by [K + ] and 
rewriting gives 



[K- 



[K- 



CxAf(7T°+7r )/ 'Ap art . 



This relation also explains that the measured ratio of K / K + is constant with 
centrality (Fig. m as the pion multiplicity does not vary with centrality. A detailed 
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Figure 10. Measured K /K + ratio, M(ir° + it )/ 'Apart and the double ratio 
([K-]/[K+])/([M(tt° + n~)]/A par t) as a function of A par t) both for Ni+Ni and 
Au+Au collisions at 1.5 A-GeV. Preliminary results. 



study for the low energies has been given in the talk by A. Forster. Figure [To| 
summarises by demonstrating the constancy of the K~/K + ratio and of the pion 
multiplicity with A part for Ni+Ni and Au+Au collisions at 1.5 A-GeV [|[ [H], |28|. 
It turns out that these ratios do not even depend on the choice of the collision 
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system. The right part of this figure exhibits the double ratio ([if~]/[X + ])/([Af (ir° + 
7r~)]/ A part ) which shows only a minor deviation from a horizontal line. This result can 
be taken as an argument that this specific channel is not far from chemical equilibrium. 
Next we test in Fig. 11 the validity of the law of mass action by plotting the 
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Figure 11. The K~/K+ ratio as a 
function of the pion multiplicity M(tt~ + 
it ) /A V art as a test of the law of mass 
action. Preliminary data. The dashed 
line shows the prediction of the statistical 
model. 



K~ /K + ratio as a function of the pion multiplicity M(ir° + Tr~)/A par t at incident 
energies from SIS up to RHIC. At SIS and AGS energies the direct relation holds, 
i.e. the K~/K + ratio rises with M(tt° + n~)/A part . At SPS and RHIC energies 
the linear relation is no longer valid; K~ are obviously produced by other channels, 
i.e. K + K~ pair production. This change of the dominating channel is well reproduced 
by the statistical model (dashed line in Fig. |TT| ) . 



2.3. Maximum relative strangeness content in heavy ion collisions around 30 A-GeV 

The experimental data from heavy ion collisions show that the K + /n + ratio rises 
from SIS up to AGS but it is larger for AGS than at the highest CERN-SPS energies 
H |l4], ^2], |3(| and decreases even further at RHIC jl5| as shown in Fig. [l^. 

This behavior is of particular interest as it could signal the appearance of 
new dynamics for strangeness production in high energy collisions. It was even 
conjectured [ ]3l| that this property could indicate an energy threshold for quark-gluon 
plasma formation in relativistic heavy ion collisions. Some transport models are able 
to describe the occurrence of this maximum within a continuous evolution of hadron 
rescattering and string degrees of freedom p^ |. 

In the following we analyze the energy dependence of strange to non-strange 
particle ratios in the framework of a hadronic statistical model. In the whole 
energy range, the hadronic yields observed in heavy ion collisions resemble those of 
a population in chemical equilibrium along a unified freeze-out curve determined by 
the condition of fixed energy /particle ~ 1 GeV || providing a relation between the 
temperature T and the baryon chemical potential fis- As the beam energy increases 
T rises and \xb is slightly reduced. Above AGS energies T exhibits only a moderate 
change and converges to its maximal value in the range of 160 to 180 MeV, while [Ib 
is strongly decreasing. 
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Figure 12. K + /ir + ratio obtained 
around midrapidity as a function of 
y/s from the various experiments. The 
dashed line shows the results of the sta- 
tistical model in complete equilibrium. 



Instead of studying the K + /ir + ratio we use the ratios of strange to non-strange 



where the quantities in angular brackets refer to the number of newly formed quark- 
antiquark pairs, i.e. it excludes all quarks that were present in the target and the 
projectile. 

Applying the statistical model to particle production in heavy ion collisions 
calls for the use of the canonical ensemble to treat the number of strange particles 
particularly for data in the energy range from SIS up to AGS j|, [54| as mentioned 
before. The calculations for Au-Au and Pb-Pb collisions are performed using a 
canonical correlation volume defined above. The quark content used in the Wroblewski 
factor is determined at the moment of chemical freeze-out, i.e. from the hadrons and 
especially, hadronic resonances, before they decay. This ratio is thus not an easily 
measurable observable unless one can reconstruct all resonances from the final-state 
particles. The results are shown in Fig. |l3| as a function of ^/s. 

The solid line (marked "sum") in Fig. [l3] describes the statistical-model 
calculations in complete equilibrium along the unified freeze-out curve |3j and with 
the energy-dependent parameters T and /is. From Fig. [l^ we conclude that around 
30 A-GeV laboratory energy the relative strangeness content in heavy ion collisions 
reaches a clear and well pronounced maximum. The Wroblewski factor decreases 
towards higher incident energies and reaches a limiting value of about 0.43. For 
details see Ref. pq| . 

The appearance of the maximum can be traced to the specific dependence of /is 
and T on the beam energy. Figure |l4| shows lines of constant \ s in the T — fj,g plane. 
As expected X s rises with increasing T for fixed fj,g . Following the chemical freeze-out 
curve, shown as a dashed line in Fig. [l4|, one can see that A s rises quickly from SIS to 
AGS energies, then reaches a maximum at [Ib ~ 500 MeV and T w 130 MeV. These 
freeze-out parameters correspond to 30 A- GeV laboratory energy. At higher incident 
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Figure 13. Contributions to the Wrob- 
lewski factor X s (for definition see text) 
from strange baryons, strange mesons, 
and mesons with hidden strangeness. 
The sum of all contributions is given by 
the full line. 
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Figure 14. Lines of constant Wrob- 
lewski factor X s (for definition see text) 
in the T — fig plane (solid lines) to- 
getherwith the freeze-out curve (dashed 
line) §. 



energies the increase in T becomes negligible but \ib keeps on decreasing and as a 
consequence A s also decreases. 

The importance of finite baryon density on the behavior of A s is demonstrated 
in Fig. [l3| showing separately the contributions to (ss) coming from strange baryons, 
from strange mesons and from hidden strangeness, i.e. from hadrons like <f> and n. 
As can be seen in Fig. ||, the origin of the maximum in the Wroblewski ratio can 
be traced to the contribution of strange baryons. This channel dominates at low y/s 
and loses importance at high incident energies. Even strange mesons exhibit a broad 
maximum. This is due to the presence of associated production of e.g. kaons together 
with hyperons. 

The energy dependence of the K + /tt + ratio measured at midrapidity is shown 
in Fig. |l2|. The model gives an excellent description of the data, showing a broad 
maximum at the same energy as the one seen in the Wroblewski factor. In general, 
statistical-model calculations should be compared with 47r-integrated results since 
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strangeness does not have to be conserved in a limited portion of phase space. A 
drop in this ratio for Air yields has been reported from preliminary results of the 
NA49 collaboration at 158 A-GeV [ p9[ . This decrease is, however, not reproduced by 
the statistical model without further modifications, e.g. by introducing an additional 
parameter j s ~ 0.7 This point might be clearer when data at other beam energies 
will become available. 

3. Summary 

Strange particle production in heavy ion collisions over a rather broad range of incident 
energies can be described by a statistical model. The production of strange particles 
close to threshold requires a canonical formulation, i.e. local strangeness conservation. 
This approach is able to explain many features of K + and K~ production at SIS 
energies. 

While for K + production it remains open whether and how chemical equlibrium 
can be reached, the situation is quite different for K~ . It is shown that the strangeness 
exchange process irA ^ N + K~ is the dominant channel for K~ production at SIS 
and likely also at AGS energies. This is demonstrated by applying the corresponding 
law of mass action. Theoretical studies confirm this interpretation. 

Using the energy dependence of the parameters T and /ib we have shown that 
the statistical-model description of relativistic heavy ion collisions predicts that the 
yields of strange to non-strange particles reaches a well defined maximum near 30 GeV 
lab energy. It is demonstrated that this maximum is due to the specific shape of the 
freeze-out curve in the T — fiB plane. In particular a very steep decrease of the baryon 
chemical potential with increasing energy causes a corresponding decline of relative 
strangeness content in systems created in heavy ion collisions above lab energies of 30 
GeV. The saturation in T, necessary for this result, might be connected to the fact 
that hadronic temperatures cannot exceed the critical temperature T c ~ 170 MeV for 
the phase transition to the QGP as found in solutions of QCD on the lattice. 

In spite of the apparent success of the statistical models it should not appear 
the impression that these models describe everything. They describe yields, particle 
ratios. Looking at spectral shapes already the expansion dynamics shows up. The 
distribution of the particles in space is a very informative quantity as e.g. |36|] . The 
description of this quantity is beyond statistical models. 

It is a pleasure for me to thank for the stimulating collaboration with P. 
Braun-Munzinger, J. Cleymans, K. Redlich, and the whole KaoS Crew (I. Bottcher, 
A. Forster, E. Grosse, P. Koczoh, B. Kohlmeyer, F. Laue, M. Menzel, L. Naumann, 
F. Piihlhofer, E. Schwab, P. Senger, Y. Shin, H. Strobele, F. Uhlig, A. Wagner, 
W. Walus). 
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